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ABSTRACT 

We study an accretion disk in which three different regions may coexist: MHD turbulent 
regions, dead zones and gravitationally unstable regions. Although the dead zones are stable, 
there is some transport due to the Reynolds stress associated with waves emitted from the 
turbulent layers. We model the transport in each of the different regions by its own a parameter, 
this being 10 to 10 3 times smaller in dead zones than in active layers. In gravitationally unstable 
regions, a is determined by the fact that the disk self-adjusts to a state of marginal stability. We 
construct steady-state models of such disks. We find that for uniform mass flow, the disk has to 
be more massive, hotter and thicker at the radii where there is a dead zone. In disks in which 
the dead zone is very massive, gravitational instabilities are present. Whether such models are 
realistic or not depends on whether hydrodynamical fluctuations driven by the turbulent layers 
can penetrate all the way inside the dead zone. This may be more easily achieved when the ratio 
of the mass of the active layer to that of the dead zone is relatively large, which in our models 
corresponds to a in the dead zone being about 10% of a in the active layers. If the disk is at 
some stage of its evolution not in steady-state, then the surface density will evolve toward the 
steady-state solution. However, if a in the dead zone is much smaller than in the active zone, 
the timescale for the parts of the disk beyond a few AU to reach steady-state may become longer 
than the disk lifetime. Steady-state disks with dead zones are a more favorable environment for 
planet formation than standard disks, since the dead zone is typically 10 times more massive 
than a corresponding turbulent zone at the same location. 

Subject headings: accretion, accretion disks — planetary systems: protoplanetary disks — stars: pre- 
main sequence 



1. Introduction 

Protoplanetary disks are believed to have a 
complex structure. In the parts which are ionized 
enough (active zones), the magnetic field couples 
to the gas and the magnetorotational instability 
(MRI; Balbus & Hawley 1991) develops, leading 
to turbulence and angular momentum transport. 
Some parts however are too cold and dense for 
the ionization to reach the level required for MHD 
turbulence to be sustained. However, in these so- 
called dead zones, there may still be some low 
level of transport due to Reynolds or non turbulent 
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Maxwell stresses (Fleming & Stone 2003, Turner 
& Sano 2008). Finally, parts of the disk may be 
dense enough that gravitational instabilities de- 
velop, redistributing mass and angular momentum 
in such a way that the disk settles into a state of 
marginal stability. To date, there is no numerical 
simulations of such global disks. Global three di- 
mensional MHD simulations so far have modeled 
stratified and fully turbulent disks, but with no 
dead zones (Fromang & Nelson 2006). 

It is of interest to study whether non-uniform 
disks composed of different regions as described 
above can be in steady-state or not. Previous 
studies (Gammie 1999, Armitage et al. 2001) 
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assumed that no transport at all took place in 
the dead zone. Therefore mass could only pile- 
up there until gravitational instabilities would de- 
velop. If the level of transport is non zero in the 
dead zone however, the disk may be able to adjust 
to reach a steady-state. It is this possibility that 
we investigate in this paper. 

In determining the characteristics of the disk 
vertical strcuture, we ignore reprocessing of the 
stellar radiation. Only heating associated with the 
different (magnetic and hydrodynamic) stresses is 
considered. Reprocessing has been shown to be 
important when the disk surface is flared (Kenyon 
& Hartmann 1987, Chiang & Goldreich 1997, 
D'Alessio et al. 1998, Dullemond et al. 2001). It 
provides an extra heating source that dominates 
the disk beyond a few AU and has important con- 
sequences on the structure of the dead zone (Mat- 
sumura & Pudritz 2006). However, in the mod- 
els we present below, the parts of the disk beyond 
1 AU are shielded from the stellar radiation by the 
inner parts. These models, which do not include 
reprocessing, are therefore self-consistent. 

The plan of the paper is as follows. In section [2] 
we describe the way we model the transport in a 
disk with turbulent regions, dead zones and gravi- 
tationally unstable regions. In section[3]we present 
the equations governing the disk vertical structure 
and describe how they are solved. In section [4] we 
present the results of the calculations, and discuss 
them in section [S] We find that steady-state so- 
lutions exist and correspond to disks which are 
thicker, hotter and more massive than disks with- 
out dead zones. 

2. Modeling of the disk transport 

Protoplanetary disks are almost certainly mag- 
netized, as they form from molecular clouds which 
are observed to contain a magnetic field. If the 
magnetic pressure is smaller than the thermal 
pressure, then the parts of the disk where the field 
couples to the gas are prone to the MRI. Its non- 
linear developement is turbulence that transports 
angular momentum outward and therefore enables 
accretion (Hawley et al. 1995, Brandenburg et al. 
1995). Recent numerical simulations of the MRI in 
shearing-boxes have shown that the level of trans- 
port depends on the Prandtl number (ratio of the 
viscosity to the resistivity; Fromang et al. 2007, 



Lesur & Longaretti 2007). The simulations, which 
have not converged yet, also show that the turbu- 
lence is not sustained when the net magnetic flux is 
zero and for the range of parameters investigated 
(relatively large Prandtl numbers). It is reason- 
able to suppose that the magnetic field in proto- 
planetary disks does not vary very much on the 
scale of a few scale heights. For this reason, and if 
these disks were to be described locally by a shear- 
ing box, it would be more appropriate to assume 
a finite net flux. In the parts of the disk that are 
coupled to the field, we will therefore assume that 
turbulence is sustained and quantify the angular 
momentum transport using the parametrization of 
Shakura & Sunyaev (1973), with olt being the ra- 
tio of the stress to the pressure (the subscript 'T' 
refers to turbulence). It has been shown that this 
prescription, which is formally valid only if the 
disk evolves under the action of a viscosity, can 
be used to describe the global properties and evo- 
lution of MHD turbulent disks, although details 
and stability issues cannot be studied within this 
framework (Balbus & Papaloizou 1999). Given 
that ax as measured in the simulations depends 
on the Prandtl number, and that Prandtl numbers 
typical of protoplanetary disks are way too small 
to be investigated numerically, we will assign val- 
ues to aT based on observations. The diffusion 
timescale is indeed limited by the disk lifetime, 
which is observed to be of a few 10 6 years. This im- 
plies «t ~ 10 -3 -10 -2 (at least in the outer parts 
of the disk which are fully turbulent). 

In protoplanetary disks, the minimum ioniza- 
tion fraction which is required for MHD turbu- 
lence to be initiated and sustained is on the order 
of 10~ 13 on scales of 1 AU (e.g., Balbus & Haw- 
ley 2000). Although very small, this level of ion- 
ization is not reached in some parts of the disks, 
which are relatively cold and dense. The extent 
of the dead zone has been computed by several 
authors, taking into account ionization by cosmic 
rays (Gammic 1996, Sano et al. 2000 who also in- 
cluded radioactivity) and X-rays (Igea & Glass- 
gold 1999, Fromang et al. 2002, Matsumura & 
Pudritz 2003 who also included cosmic rays and 
radioactivity). In this paper, we will denote by £t 
the surface density of the active layer, i.e. the sur- 
face density beyond which ionizing photons cannot 
penetrate. Depending on the energy of these pho- 
tons, St may vary in the range 10-100 g cm -2 
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(Gammie 1996, Fromang et al. 2002). Note that 
the ionization in the disk inner parts is thermal, 
with the most abundant ions being Na + and K+. 
For the ionization fraction to be larger than 10~ 13 , 
the temperature has to be higher than some value 
Ti which is close to 10 3 K over a large range of 
densities (Balbus & Hawley 2000). There is evi- 
dence from numerical simulations that some level 
of transport can be present in a dead zone that 
is sandwiched in between turbulent layers (Flem- 
ing & Stone 2003), as the turbulence drives hydro- 
dynamic waves that propagate in the dead zone 
and are associated with a Reynolds stress. More 
recently, Turner & Sano (2008) have also proposed 
that magnetic field be able to diffuse from the ac- 
tive layers into the dead zone. Although the re- 
sistivity there is too high for MHD turbulence to 
be sustained, there is a large scale magnetic stress 
associated with the field that drives accretion. (If 
dust grains are present in the disk however, the 
magnetic field may not able to diffuse down to the 
disk midplane at all radii.) The transport and 
energy deposition in the dead zone due to these 
stresses are non local. However, for illustrative 
purposes, we will parameterize this transport with 
a Shakura & Sunyaev parameter that we will de- 
note an (where the subscript 'D' refers to dead 
zone). Fleming & Stone (2003) found that the 
Reynolds stress in the dead zone is about 10% of 
the Mawell stress in the active layers. Normalized 
by the thermal pressure, this translates into 
being between 0.07 and 0.3 times olt- However, 
in their simulations, the mass of the dead zone 
was only a few times that of the active layer. We 
might expect lower values of the Reynolds stress 
deep into the dead zone when the mass of the ac- 
tive layer is significantly smaller than that of the 
active zone, which is the case in the models we 
present below. Therefore, we will take an in the 
range 10 _3 o:t - 10 _1 aT- 

In the parts of the disk which are dense and 
cold, gravitational instabilities develop. Spiral 
density waves grow and transport angular mo- 
mentum outward and mass inward. The level 
of saturation is determined by the fact that the 
disk settles into a state of marginal stability, with 
the Toomre parameter Q being on the order of 
1.5 (Laughlin & Bodcnheimcr 1994, Lodato & 
Rice 2004). It was shown by Balbus & Pa- 
paloizou (1999) that in general, in a gravitation- 



ally unstable disk, the transport is non local and 
cannot be described using the viscous disk theory. 
However, these authors pointed out that a local 
transport model might apply if the disk maintains 
itself in a state of marginal stability, which is the 
case considered here. In the regions where the 
disk is gravitationally unstable, we will therefore 
use the Shakura & Sunyaev prescription and will 
adjust the transport parameter, that we will de- 
note «g (where the subscript 'G' refers to gravity), 
such as to have Q ~ 1.5. 

It is important to keep in mind that here the 
disk is modeled in a very schematic way. We 
believe this does not affect the gross features of 
the models we present below, although the details 
should be considered with caution. 

3. Disk vertical structure 

We consider a system of cylindrical coordinates 
(r, ip, z) based on the central star, with z = being 
the disk midplane. P is the pressure, p is the mass 
density per unit volume, T is the temperature, Q is 
the angular velocity, v is the 'enhanced' kinematic 
viscosity associated with ax in turbulent layers, 
afl in dead zones and <xq in gravitationally unsta- 
ble regions, and k is the opacity, which in general 
depends on both p and T. The thin disk approxi- 
mation is used throughout (and checked to be valid 
a posteriori), so that U 2 = GM±/r 3 , M± being the 
mass of the central star and G the gravitational 
constant. 

The disk vertical structure is described by the 
equation of vertical hydrostatic equilibrium: 

p oz 

together with the energy equation, which states 
that the rate of energy removal by radiation is 
locally balanced by the rate of energy production 
by viscous dissipation: 

BF 9 9 

°jL = \ pv tf = » Q nP, (2) 

where F is the radiative flux of energy through 
a surface of constant z. To write the last term 
of this equation, we have used v = ac 2 /Sl, where 
c s = (P/p) 1 / 2 is the isothermal sound speed and 
a = a T in active layers, a = a D in dead zones 
and a — a G in the gravitationally unstable parts 
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of the disk. Therefore a is a function of both z 
and r. The flux F is given by: 



F = 



-16ctT 3 dT 
Sup dz ' 



(3) 



with a being the Stefan-Boltzmann constant. In 
principle, equation ^ is valid only at radii where 
the disk is optically thick. However, when the disk 
is optically thin, i.e. when np integrated over the 
disk thickness is small compared to unity, the tem- 
perature gradient given by equation ((3|) is small, 
so that the results we get are consistent in that 
case also. To close the system of equations, we 
adopt the equation of state of an ideal gas: 



P = 



pkT 
2m H ' 



(4) 



where k is the Boltzmann constant and 2m# is 
the mass of the hydrogen molecule, which is the 
main component of protostellar disks at the tem- 
peratures and densities of interest here. 

3.1. Boundary conditions 

We have to solve three first order ordinary dif- 
ferential equations for the three variables F, P (or 
equivalently p), and T as a function of z at a given 
radius r. Accordingly, we need three boundary 
conditions at each r. These have been described 
in detail in Papaloizou & Terquem (1999), so here 
we just recall briefly their expression. We denote 
with a subscript s values at the disk surface. The 
flux at the surface is given by: 



f s = —Km 2 , 

8tt 



(•5) 



where M = Ztt{v)Yj, with £ = J^ H pdz being the 

disk surface mass density and (y) = J H R pvdz/Y, 
being the vertically averaged viscosity. If the disk 
is in a steady state, M does not vary with r and 
is the constant accretion rate through the disk. 
H and —H are the upper and lower boundaries 
of the disk, respectively. The surface pressure is 
given by: 



P, 



Q 2 Ht, 



a b 



(6) 



where T a f, is the optical depth above the disk. 
Since we have defined the disk surface such that 



the atmosphere above the disk is isothermal, we 
have to take T a & <C 1. Providing this is satis- 
fied, the results do not depend on the value of 
T a b we choose (see Papaloizou & Terquem 1999). 
Finally, the surface temperature satisfies the fol- 
lowing equation: 



2a (Tf - Tu) 



9a s kT s Cl 
8/im_f//« s 



8tt 



-MCI 



0. (7) 



Here the disk is assumed immersed in a medium 
with background temperature Tf,, so that the sur- 
face temperature remains finite. 

3.2. Model Calculations 

We look for steady-state solutions, so that M 
does not depend on r. At a given r we then solve 
equations flU , © and fl3J) with the boundary con- 
ditions ©, © and to find the dependence of 
the state variables on z. The (Rosseland) opac- 
ity is taken from Bell & Lin (1994). This has 
contributions from dust grains, molecules, atoms 
and ions. It is written in the form k = Kip a T b 
where m, a and b vary with temperature. To 
carry out the integration, we start from an esti- 
mated value of H and iterate until the condition 
F = at z — is satisfied (see details in Pa- 
paloizou & Terquem 1999). The integration starts 
at z = H, where we have a — ay, assuming the 
disk is not gravitationally unstable, and proceeds 
all the way down to z = 0. At each value of z, 
we compute £(z) = / p{z')dz' and compare it 
with S T . When S(z) > S T , and if T(z) < T, 
we fix a = Old- Once the vertical integration 
is finished, we calculate the Toomre parameter 
Q = c s CI/ttGT,, where S is now the total sur- 
face density and c s is evaluated at the midplane. 
When Q < 1.5, we increase a, that we now call etc, 
slightly, assume that it is now independent of z, 
and redo the integration. We iterate the calcula- 
tion, increasing ac at each step, until Q becomes 
larger than 1.5. As the step with which we in- 
crease aa has a finite size, we usually end up with 
a value of Q between 1.5 and 2 after at most 20 
iterations. 

In the calculations presented here, we have 
taken M* = 1 M , T = 10 3 K, the optical depth 
of the atmosphere above the disk surface r a f, = 
10 -2 and a background temperature Tf, — 10 K. 
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In the optically thick regions of the disk, the value 
of H is independent of the value of r a b we choose. 
This is not the case in optically thin regions where 
we find that, as expected, the smaller T a &, the 
larger H. However, this dependence of H on r a f, 
has no physical significance, since the surface mass 
density, the optical thickness through the disk and 
the midplane temperature hardly vary with r a b- 
This is because the mass is concentrated towards 
the disk midplane in a layer with thickness inde- 
pendent Of Tab- 

4. Results 

In figure [TJ we plot H/r versus r (from 0.05 
to 100 AU) for M = 10~ 8 M© yr" 1 , a T = 10~ 2 , 
a D = 10~ 4 and E T = 10 and 50 g cm~ 2 . The 
solid line delimits the surface of the disk whereas 
the dotted lines delimit the different regions in the 
disk. For Ex = 10 g cm~ 2 , there is a dead zone 
that extends from about 0.1 AU to about 6 AU. 
Note that at small radii, this dead zone is present 
only at some intermediate values of z, although 
the ionizing photons cannot penetrate down to 
the disk midplane. This is because, if the disk 
starts turbulent there, the temperature near the 
midplane is large enough for the gas to be ther- 
mally ionized and therefore to remain turbulent. 
And if the disk does not start turbulent there, the 
temperature in the dead zone is large enough for 
thermal ionization to be above the citical level for 
the onset of the MM, so that any perturbation will 
lead to turbulence. From r ~ 6 AU up to about 
27 AU, the disk is gravitationally unstable. In this 
region, «g increases with r from 10~ 4 up to 10 -2 . 
The disk is indeed cooler in the outer parts so 
that, if o>g were kept constant, Q would decrease 
when r increases. Therefore, to maintain Q con- 
stant, Qg has to increase with radius. The inner 
parts of the disk, where the temperature is larger 
than Ti, and the surface and outer parts, where 
the ionizing photons can penetrate, are turbulent. 
For St = 50 g cm -2 , the dead zone extends from 
r ~ 0.2 AU to r ~ 5 AU and there is no grav- 
itationally unstable region. This is because the 
dead zone is less extended and less massive than 
for smaller values of Ey. As the ionizing photons 
can penetrate deeper than in the previous case, 
the turbulent zone is much more extended. 

In this model, there is a range of radii where 



the mass of the dead zone is much larger than 
that of the active layer. We denote E a and E^ 
the column density of the active layer and dead 
zone, respectively. At r = 1 AU, where the dead 
zone is vertically very extended, we have E a /Ed = 
0.6 xl0~ 2 and4xl0~ 2 for E T = 10 and 50 gem- 2 , 
respectively. Whether hydrodynamic fluctuations 
driven by the turbulence in the active layers can 
penetrate down to the disk midplane in that case 
is not clear. 

In order to get a model with higher values of 
£ a /£d, we have investigated the case ax = 10 -2 
and a D = 10~ 3 (M = 10 -8 M Q yr" 1 being the 
same as above). This is close to the values found 
by Fleming & Stone (2003) in their simulations 
in which E a /Ed was varied between 0.187 and 
0.75. Figure [5] is the same as figure [1] but with 
this larger value of old- We note that the range 
of radii at which the disk is not fully turbulent 
is the same as before. Here however, there is no 
gravitationally unstable region. The dead zone be- 
ing less massive, it remains gravitationally stable. 
For Yjt = 10 g cm~ 2 , the dead zone now extends 
from r - 0.1 AU up to about 27 AU. The ver- 
tical extent of the dead zone is smaller than in 
the case cud — 10~ 4 . At r = 0.4 AU, where the 
dead zone is the more extended, E a /Ed = 10~ 2 
and 0.1 for = 10 and 50 g cm -2 , respectively. 
This latter value is closer to those investigated by 
Fleming & Stone (2003). Note that if we were to 
assume that the Reynolds stress in the midplane 
were 10% of the Maxwell stress at the boundary 
of the dead/active zone, as found by Fleming & 
Stone (2003), we would get a_o = 0.03o;t after 
normalizing the stress by the thermal pressure at 
these locations. But of course the angular momen- 
tum flux is spread throughout the different zones 
and in principle a should vary with z. 

In figure 02 we plot H/r, E and the midplane 
temperature T m versus r for the same parameters 
as in figures Q] and O For comparison, we also plot 
the curves corresponding to standard disk models 
(with no dead zone and no adjustment in the grav- 
itationally unstable parts) with constant a = cut 
and a = old. As expected, in the parts where 
there is a dead zone, the disk is thicker, more mas- 
sive and hotter. This is because the vertically av- 
eraged value of a decreases as the vertical extent 
of the dead zone increases, and H/r, E and T m 
increase as a decreases (e.g., Frank et a!. 1992). 



5 



We see that for low values of St, i-e. when the 
ionizing photons do not penetrate deep into the 
disk, the turbulent layer is thin and the proper- 
ties of the disk at the radii where there is a dead 
zone are very close to that of a standard disk with 
constant ol = old- For all the values of St, H/r 
has a maximum at a fraction of or around 1 AU, 
which produces a puffing up of the disk. This is 
because the mass density is very high there, so the 
dead zone is vertically extended, which makes the 
disk hotter. Note that a maximum is also present 
in the standard disk models with constant a. The 
profiles of H/r, S and T m tend to be a bit steeper 
than in standard disk models with no dead zone, 
as they are bracketed by the profiles corresponding 
to constant a — or and a = old models. 

In figure [H we plot p and T as a function of z/r 
at r = 1 AU for a D = 1(T 4 and at r = 0.4 AU 
for Od = 1(U 3 for the same parameters as in fig- 
ure [3] When the dead zone extends down to the 
disk midplane, the disk structure is very close to 
that of a standard disk with constant a = ao for 
the lowest value of St- For the largest values of 
St, the dead zone is small and the disk is simi- 
lar to a standard disk with constant a = olt- For 
old — 1CP 4 , we actually cannot distinguish be- 
tween the two models for St = 10 2 g cm~ 2 . We 
see that the mass density in the dead zone is typi- 
cally 10 times larger than that of a standard fully 
turbulent disk with constant a — or- 

In figure we plot H/r versus r for M = 
10~ 8 M Q yr" 1 , a T = W 3 , a D = 10~ 5 and 
St = 50 g cm~ 2 , i.e. the same values as above 
except for lower values of or and a d ■ We checked 
that decreasing or and old by a factor of ten has 
almost exactly the same effect as increasing M by 
the same factor, all the other parameters being 
kept fixed. When or and old are smaller or M is 
larger, the disk is more massive and therefore grav- 
itational instabilities develop in a more extended 
region, up to the disk outer boundary. To keep Q 
constant ~ 1.6 there, etc has to reach 2.5 x 10~ 3 
in the disk outer parts, i.e. to become larger than 
or- 

We have also run models with ar = 10 -2 and 
old = 10 _3 «t- We get similar results in that case, 
with the locations where there is a dead zone being 
hotter, thicker and more massive. 

We note that in all the models presented above, 
the puffing up of the inner parts of the disk casts 



a shadow over the outer parts, which therefore do 
not reprocess stellar radiation. The models pre- 
sented here are self-consistent in a sense that if 
the disk is not illuminated by the star to begin 
with, it will not be able to reprocess the stellar 
radiation at any time. There may however exist 
another type of solutions in which stellar irradi- 
ation is important. Indeed, if the outer parts of 
the disk are illuminated at some point (before a 
steady-state is reached for instance), a flaring is 
produced that could be maintained. To test this 
hypothesis, we have modified the surface temper- 
ature. Instead of using equation to calculate 
T s , we fix T s = 300 (r/AU) _1/2 K. This corre- 
sponds to T s ~ 950 K, 300 K and ~ 95 K at 
r = 0.1, 1 and 10 AU, respectively. These values 
are close to the surface temperatures computed by 
D'Alessio et al. (1999) in irradiated disk models 
for M = 10~ 8 Mq y^ 1 and a = 10~ 3 , this value 
of a being intermediate between «t and aD ■ With 
this new T s , the profile of H/r we get still shows a 
shadowing of the outer parts. This suggests that 
even if the disk outer parts are illuminated by the 
central star at some point, before a steady-state is 
reached, they cannot flare enough to remain irra- 
diated after a steady-state is established, so that 
reprocessing of stellar radiation does not play a 
role in the steady-state models. 

5. Discussion and conclusion 

In this paper, we have considered a disk in 
which transport is produced by either: i) Maxwell 
stress resulting from MHD turbulence in the (grav- 
itationally stable) regions where the gas couples 
well to the magnetic field; ii) gravitational stress 
in the parts which are gravitationally unstable; 
iii) Reynolds stress associated with hydromagnetic 
waves driven by the turbulence in adjacent lay- 
ers and Maxwell stress associated with large scale 
field siphoned off from these turbulent layers in 
the parts which are not turbulent and are gravi- 
tationally stable (dead zones). We have modeled 
the transport using the Shakura & Sunyaev pre- 
scription, with old being 10 to 10 3 times smaller 
than olt, and olq being adjusted such as to give 
a Toomre parameter Q ~ 1.5 in the gravitation- 
ally unstable regions. Although such a modeling is 
in principle not valid when the transport is non- 
local, we believe it does not invalidate the gross 
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features of the models we have presented. 

We have found that steady-state models of such 
a disk exist, and that they are physically reason- 
able. Since M = 37r(^)E oc ac s HY,, mass flow 
through the disk can be uniform only if the disk 
is more massive, hotter and thicker at the radii 
where the vertically-averaged value of a is smaller, 
i.e. where there is a dead zone. Note that even 
though the disk is thicker at the locations where 
there is a dead zone, it remains thin for the param- 
eters investigated here. In disks in which the dead 
zone tends to be massive, which is the case for 
the lowest values of oid and/or the lowest values 
of £t investigated here, gravitational instabilities 
control the dynamics of part of the disk. 

Whether these models are realistic or not de- 
pends on whether hydrodynamical fluctuations 
driven by the turbulent layers can penetrate all 
the way inside the dead zone. This may be more 
easily achieved when the ratio of the mass of the 
dead zone to that of the active layer is the smallest, 
which in our models corresponds to old j cut = 0.1. 

If the disk is at some stage of its evolution out of 
steady-state, then the surface density will change 
in such a way as for the disk to evolve toward 
steady-state. To see this, let us consider a disk in 
which M decreases inward at some location. Then 
mass will accumulate there as accretion through 
this region is slower, until S is large enough for 
the flow to be steady. In contrast, if M increases 
inward at some location, accretion there will be 
faster and mass will be depleted until £ is re- 
duced to the level where a steady-state is reached. 
The timescale t v for establishing steady-state at 
the radii where there is a dead zone may be long 
though, as it is given by: 



where (a) — f_ H padzfE. With (a) = 10 2 and 
H/r = 0.1, we get t„~5x 10 2 and 6 x 10 3 years 
at 1 and 5 AU, respectively, which is much smaller 
than the disk lifetime. However, if the dead zone 
at these radii is vertically very extended, we have 
(a) ~ old , so that the viscous timescale can get 
100 times longer if a D = 10~ 4 . At 1 AU we 
still get a timescale significantly smaller than the 
disk lifetime, but this is not the case at 5 AU, 
where the disk may not be able to reach steady- 



state. If these parts of the disk build-up from 
mass inflowing from further out, then the mass 
in the dead zone will slowly increase. This pro- 
cess will never produce outbursts as envisioned by 
Gammie (1999) though. Indeed, before enough 
mass could pile-up for an outburst to be produced, 
a steady-state would be achieved in which mass 
would be transported either by density waves or 
Reynolds stress driven by the turbulent layers. 

Note that the spectral energy distribution of a 
steady disk such as those considered in this paper 
is not different from that of a standard disk with 
no dead zone, as the total flux emitted at the sur- 
face depends only on M. The surface temperature 
also does not depend on Sr, as can be seen from 
equation ([7]). However, the temperature just be- 
low the surface is significantly larger in disks with 
dead zones than in standard disks. The chemistry 
there would therefore be different. Grain prop- 
erties may also differ. To study this, irradiation 
of the disk by the central star would have to be 
taken into account. Although it will not be im- 
portant beyond the regions where there is a dead 
zone, as these parts are shielded from the stellar 
radiation due to the puffing up of the disk, it will 
be important in the parts where H/r increases. 

Finally, we comment that a steady-state disk 
with a dead zone is a more favorable environment 
for planet formation than a standard disk, as the 
dead zone, which in general encompasses the re- 
gion of planet formation, is significantly more mas- 
sive. Planet formation would therefore be faster 
there (Lissauer 1993). Note however that there is 
an issue as how to stop type I migration of the 
cores in these dead zones. Turbulence, which has 
been shown to alter this type of migration (Nel- 
son & Papaloizou 2004), cannot be invoked. But 
if large scale fields are siphoned off in these re- 
gions from the turbulent layers, they may pre- 
vent the cores from migrating too fast (Terquem 
2003, Fromang et al. 2005, Muto et al. 2008). 
Planetary cores may also be stopped at the inter- 
face between the dead zone and the inner turbu- 
lent region, where the surface mass density varies 
sharply, as suggested by Masset et al. (2006). 
This process however relies on the (positive) coro- 
tation torque becoming dominant over the Lind- 
blad torque due to the density gradient and was 
studied in the context of a laminar viscous disk. It 
is not clear how the corotation torque is affected 
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by the turbulence at the transition between the 
active and dead zones, and whether it would still 
be as efficient in this context. 
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cal Physics at UCSB, funded by the US National 
Science Foundation under grant PHY05-51164. I 
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Fig. 1.— H/r versus r (in au) for M = 10" 8 M Q yr" 1 , a T = 10" 2 , a D = 10" 4 , Tj = 10 3 K and 
St = 10 g cm~ 2 (upper plot) and 50 g cm~ 2 (lower plot). The solid line delimits the surface of the disk 
whereas the dotted lines delimit the different regions in the disk. For St = 10 g cm~ 2 , there is a dead 
zone that extends from r ~ 0.1 au to r ~ 6 au and the disk is gravitationally unstable from r ~ 6 au up 
to r ~ 27 au. For Sy = 50 g cm~ 2 , the dead zone extends from r ~ 0.2 au to r ~ 5 au and there is no 
gravitationally unstable region. 
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M=10- 8 Mgyr- 1 , a T = 1 0~ 2 , a D = 1 0~ 3 





Fig. 2. — Same as Fig. Q] but for ar = 10 -2 and an = 1CP 3 . For = 10 g cm -2 , the dead zone extends 
from r ~ 0.1 au to r ~ 27 au whereas for St = 50 g cm -2 , it extends from r ~ 0.2 au to r ~ 5 au. There is 
no gravitationally unstable region in that case and the vertical extent of the dead zone is smaller than for 
a D = 10" 4 . 
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Fig. 3. — H/r (upper plots), S in g cm~ 2 (middle plots) and midplane temperature T m in K (lower plots) 
versus r (in au) for M = 10~ 8 M yr _1 , olt = 10~ 2 and q.d = 1CT 4 (left plots) and cut = 10~ 2 and 
olb = 1CP 3 (right plots). The solid and dotted-dashed curves correspond to standard disk models (no dead 
zones) with constant a = ax and a = an, respectively. The dotted, long-dashed and short-dashed curves 
correspond to the disk models with St = 10 2 , 50 and 10 g cm~ 2 , respectively. In the parts where there is a 
dead zone, the disk is thicker, more massive and hotter. 
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Fig. 4. — p in g cm -3 (upper plots) and T in K (lower plots) as a function of z/r for M = 10~ 8 yr -1 , 
or = 10~ 2 and old = 10~ 4 at r = 1 au (Ze/t plots) and or = 10~ 2 and od = 10~ 3 at r = 0.4 au (rig/?i plots). 
The solid and dotted-dashed curves correspond to standard disk models (no dead zones) with constant 
a = a? and a ~ au, respectively. The dotted, long-dashed and short-dashed curves correspond to the disk 
models with St = 10 2 , 50 and 10 g cm~ 2 , respectively. On the left plots, the solid and dotted lines are 
indistinguishable. The vertical structure of the disk depends on the vertical extent of the dead zone. 
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M=10- 8 Mgyr- 1 , a T = 1 0~ 3 , a D = 1 0~ 5 




Fig. 5. — Same as Fig.Q]but for ar — 10 3 , old = 10 5 and E = 50 g cm 2 . Here the disk is gravitationally 
unstable from r~3au all the way up to the outer radius. 
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